Methods for Differential Image Quality Enhancement for a Multiple Detector System, Systems and Use Thereof

ABSTRACT

Methods for differential image quality enhancement for a detection system including multiple electromagnetic radiation detectors which include obtaining an image from a chemical band electromagnetic radiation detector and an image from a reference band electromagnetic radiation detector. Each of the images includes a plurality of pixels, each pixel having an associated intensity value. One or more intensity values of a plurality of pixels from the reference band image are adjusted based on one or more intensity value parameters of the chemical band image.

CROSS-REFERENCE TO RELATED APPLICATION

This application claims the benefit of U.S. Provisional Application No. 62/015,921, filed Jun. 23, 2014, entitled METHODS FOR DIFFERENTIAL IMAGE QUALITY ENHANCEMENT FOR A MULTIPLE DETECTOR SYSTEM, SYSTEMS AND USE THEREOF, the entirety of which is incorporated by reference herein.

FIELD

The present disclosure relates to systems and methods for differential image quality enhancement of images obtained from multiple detector systems and the use thereof. In particular, the multiple detector systems may be used to identify chemical species, for example detecting chemical emissions.

BACKGROUND

This section is intended to introduce various aspects of the art, which may be associated with one or more embodiments of the present disclosure. This discussion is believed to assist in providing a framework to facilitate a better understanding of particular aspects of the present disclosure. Accordingly, it should be understood that this section should be read in this light, and not necessarily as admissions of prior art.

Chemical usage is a fundamental aspect of current civilization. Facilities for the production, processing, transportation, and use of chemical species continue to be built in locations around the world. Thus, detection of chemical species is a continuing focus.

An example of chemical species detection is gaseous leak detection. As the efficiency of facilities becomes increasingly important, even minor losses of chemical species such as hydrocarbons can add to cost or create issues for regulatory agencies.

Hydrocarbons in facilities may be lost due to process limitations or process upsets leading to flaring or leaks. While some of these issues can be directly improved by design, leaks still provide a challenge, as they may occur on any number of different process equipment types. For example, leaks can originate from pipe flanges, valves, valve stems, sampling systems, and any number of other locations. As equipment is used and ages, leaks become increasing probable.

Conditions within a facility can increase the probability of leakage or exacerbate leaks when they do form. For example, facilities using high pressures or cryogenic temperatures can increase the probability of leaks. LNG plants are an example of such facility conditions. The number of LNG liquefaction plants around the world is growing rapidly and as these plants age, there is the increasing potential for hydrocarbon leaks to develop.

Early detection and repair of leaks can be useful in preventing any number of issues, such as increased costs and regulatory issues. Leaks may be detected by operators, for example, by visually observing the release, smelling the hydrocarbons, or hearing noise caused by the release. However, most hydrocarbon vapors are not visible to the naked eye. Further, there is often a high level of equipment congestion in plants, which may place a leak point behind another piece of equipment. In addition, hydrocarbons may have a minimal odor and, thus, may not be detected by smell. Detecting a small leak by sound is more improbable, as the very high level of ambient noise and safety equipment such as earplugs makes it unlikely that the leak will be heard.

Leak detection systems have been installed in many facilities. One such system may include combustible gas detectors that monitor the concentration or lower explosive limit (LEL) of hydrocarbon vapors at a particular location, providing a measurement of a hydrocarbon level at a point in an area. An array of point measurement systems may then be used to track a vapor release across the area. However, point detection systems may not detect small releases, such as from small leaks or new leaks, the amount of hydrocarbons released, and the like.

Another leak detection system that has been used utilizes a detection method that detects hydrocarbons in a line across a plant environment, for example, by directing a light source at one edge of an area towards a spectroscopic detector at another edge of the area. While such systems may be useful for monitoring compliance for regulatory issues, they do not necessarily identify a location of a release along the line. Further, they may not detect small releases at all for the same reasons as the point detectors, e.g., the hydrocarbons may be too dilute to detect or may be blown away from the detection line by the wind.

Another leak detection system has been described that can detect releases by imaging areas using cameras which can directly show an image of a hydrocarbon plume. One such system is described in Hackwell, J. A., et al., “LWIR/MWIR Hyperspectral Sensor for Airborne and Ground-based Remote Sensing,” Proceedings of the SPIE, Imaging Spectroscopy II, M. R. Descour, and J. M. Mooney, Eds., Vol. 2819, pp. 102-107 (1996). The system was named a spatially-enhanced broadband array spectrograph system (SEBASS). The SEBASS system was intended to explore the utility of hyperspectral infrared sensors for remotely identifying solids, liquids, and gases in a 2 to 14 micrometer spectral region often used to provide a chemical fingerprint. The SEBASS system allows the imaging and identification of chemical materials, such as hydrocarbon plumes, in an environment.

In a presentation entitled “The Third Generation LDAR (LDAR3) Lower Fugitive Emissions at a Lower Cost” (presented at the 2006 Environmental Conference of the National Petrochemical & Refiners Association, Sep. 18-19, 2006), Zeng, et al., discloses an autonomous system for leak detection that uses a camera to identify leaks in a particular area of a plant. Infrared (IR) video images from the camera are processed using software to minimize background and noise interference and the likely volatile organic compound (VOC) plumes are isolated using an algorithm. A plume index (PI) is calculated based on the number and intensity of pixels in the processed VOC plume image. If the PI is greater than an experimentally determined threshold value, an action can be triggered, such as an alarm or a video capture, for confirmation.

Another such system is described in WO2012/134796. The apparatus described therein includes multiple detectors configured to address complex interferences, such as moving equipment, people, vehicles, or steam, which can lead to false detections with a single detector system.

While the existing systems attempt to minimize background and noise interference, there is still a desire to obtain improved images for more accurate detection of chemical species.

SUMMARY

This summary is meant to provide an introduction of the various embodiments further described herein and is not meant to limit the scope of claimed subject matter.

The present disclosure relates to a method for differential image quality enhancement for a detection system including multiple electromagnetic (EM) radiation detectors. The method includes obtaining at least a first image from a first electromagnetic radiation detector of the detection system. The image includes a first plurality of pixels, each pixel having an associated intensity value. At least a second image from a second electromagnetic radiation detector of the detection system is obtained. The image includes a second plurality of pixels, each pixel having an associated intensity value. One or more intensity values of the second plurality of pixels of the second image from the second electromagnetic radiation detector are adjusted based on one or more intensity value parameters of the first image from the first electromagnetic radiation detector.

In another aspect, the present disclosure relates to a system for detecting at least one chemical species comprising: a lens; a beam splitter; a first electromagnetic radiation detector, a second electromagnetic radiation detector; and an analysis system. The analysis system includes code within a processor, a non-transitory, computer-readable medium, or a combination thereof. The code is configured to direct the processor to: identify at least a first image from the first electromagnetic radiation detector, the first image including a first plurality of pixels, each pixel having an associated intensity value; identify at least a second image from the second electromagnetic radiation detector, the second image including a second plurality of pixels, each pixel having an associated intensity value; and adjust one or more intensity values of the second plurality of pixels of the second image from the second electromagnetic radiation detector based on one or more intensity value parameters of the first image from the first electromagnetic radiation detector.

In yet another aspect, the present disclosure relates to a method of detecting at least one chemical species. The method including providing a detection system as described herein. A beam of electromagnetic radiation is split into at least a first electromagnetic radiation beam and a second electromagnetic radiation beam. At least a portion of the first electromagnetic radiation beam is received by a first electromagnetic radiation detector and at least a portion of the second electromagnetic radiation beam is received by a second electromagnetic radiation detector. At least a first image of the first electromagnetic radiation beam is generated by the first electromagnetic radiation detector and at least a second image of the second electromagnetic radiation beam is generated by the second electromagnetic radiation detector. The image quality is enhanced in accordance with the methods described herein. The presence of at least one chemical species is determined based, at least in part, on a comparison of the enhanced first image from the first electromagnetic radiation detector and the enhanced second image from the second electromagnetic radiation detector.

Other aspects of the present disclosure will be apparent from the following description and the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing and other advantages of the present disclosure may become apparent upon reading the following detailed description and upon reference to the drawings, in which:

FIG. 1 illustrates a detection system in accordance with one or more embodiments of the present disclosure.

FIGS. 2A and 2B illustrate histograms of pixel intensities of images from a multiple detector system according to one or more embodiments of the present disclosure.

FIG. 3 illustrates a flow chart for a method of enhancing image quality of a detection system including multiple EM radiation detectors in accordance with one or more embodiments of the present disclosure.

FIG. 4 illustrates an image from a chemical band EM radiation detector as well as an image from a reference band EM radiation detector and a resultant image in accordance with one or more embodiments of the present disclosure.

FIG. 5 illustrates a flow chart for a method of detecting a chemical species in accordance with one or more embodiments of the present disclosure.

DETAILED DESCRIPTION

In the following detailed description section, the specific embodiments of the present disclosure are described in connection with one or more embodiments. However, to the extent that the following description is specific to a particular embodiment or a particular use of the present disclosure, this is intended to be for exemplary purposes only and simply provides a description of the one or more embodiments. Accordingly, the disclosure is not limited to the specific embodiments described below, but rather, it includes all alternatives, modifications, and equivalents falling within the true spirit and scope of the appended claims.

Various terms as used herein are defined below. To the extent a term used in a claim is not defined below, it should be given the broadest definition persons in the pertinent art have given that term as reflected in at least one printed publication or issued patent.

Certain terms are used throughout the following description and claims to refer to particular features or components. As one skilled in the art would appreciate, different persons may refer to the same feature or component by different names. This document does not intend to distinguish between components or features that differ in name only. The drawing figures are not necessarily to scale. Certain features and components herein may be shown exaggerated in scale or in schematic form and some details of conventional elements may not be shown in the interest of clarity and conciseness. When referring to the figures described herein, the same reference numerals may be referenced in multiple figures for the sake of simplicity. In the following description and in the claims, the terms “including” and “comprising” are used in an open-ended fashion, and thus, should be interpreted to mean “including, but not limited to.”

As used herein, a plurality of items, structural elements, compositional elements, and/or materials may be presented in a common list for convenience. However, these lists should be construed as though each member of the list is individually identified as a separate and unique member. Thus, no individual member of such list should be construed as a de facto equivalent of any other member of the same list solely based on their presentation in a common group without indications to the contrary.

Concentrations, quantities, amounts, and other numerical data may be presented herein in a range format. It is to be understood that such range format is used merely for convenience and brevity and should be interpreted flexibly to include not only the numerical values explicitly recited as the limits of the range, but also to include all the individual numerical values or sub-ranges encompassed within that range as if each numerical value and sub-range is explicitly recited. For example, a numerical range of 1 to 4.5 should be interpreted to include not only the explicitly recited limits of 1 to 4.5, but also to include individual numerals such as 2, 3, 4, and sub-ranges such as 1 to 3, 2 to 4, etc. The same principle applies to ranges reciting only one numerical value, such as “at most 4.5”, which should be interpreted to include all of the above-recited values and ranges. Further, such an interpretation should apply regardless of the breadth of the range or the characteristic being described.

The term “electromagnetic radiation” or “EM radiation” shall mean electromagnetic waves or photons that carry energy from a source. EM radiation is often categorized into spectral ranges by its interaction with matter, for example radio waves, microwaves, infrared, visible light, ultraviolet light, x-rays, and gamma rays. As used herein, x-rays include wavelengths in the range of from 0.01 nanometers (nm) to 10 nm. Ultraviolet (UV) light, or the UV spectrum, includes light having wavelengths of 190 nm to 400 nm. Visible light or the visible spectrum includes light that is detectable by a human eye, for example from 400 nm to 700 nm. In the UV and visible spectral ranges, chemical species may absorb energy through electronic transitions in which an electron is promoted from a lower orbital to a higher orbital. IR light, or the IR spectrum, includes light at wavelengths longer than the visible spectrum but generally lower than the microwave region.

For example, the IR spectrum may include light having a wavelength between 0.7 micrometers (microns) and 14 microns (700 nm to 14000 nm) in length. At the longer wavelength end of this range, at 10 microns to 14 microns (the far IR), chemical species may absorb energy through rotational transitions. At an intermediate wavelength range of 2.5 microns to 10 microns (the mid-IR), chemical species may absorb energy through vibrational transitions. At the lower end of the wavelength range, at 0.7 microns to 2.5 microns (the near-IR), chemical species may absorb energy through vibrational transitions and through similar processes as visible and UV light, e.g., through electronic transitions. EM radiation detectors form images from EM radiation in the visible spectrum, IR spectrum, or UV spectrum.

The term “camera” as used herein means a device that can obtain a two dimensional image or a sequence of two dimensional images or frames (such as video or series of still images) in a particular EM radiation spectral range.

The term “chemical species” is any compound that may be released into the environment as a gas, vapor, or liquid. Examples of chemical species that may be detected using the systems and methods described herein include both hydrocarbons and other chemical species. Chemical species that may be detected include but are not limited to hydrocarbon gas or vapors released in a cloud or plume into the air at an LNG plant or other facility or hydrocarbon liquids (e.g., oil) forming a slick on top of a body of water. Non-hydrocarbon species that may be detected include but are not limited to hydrogen fluoride gas released in a refinery, chlorine gas released in a water treatment facility, or any number of other liquids, gases, or vapors. A chemical species may also be deliberately added to a process stream to enhance the detection of a plume using the systems and methods described herein.

The term “facility” as used herein means any location including a tangible piece of physical equipment. For example, a tangible piece of physical equipment through which hydrocarbon fluids are produced from a reservoir, injected into a reservoir, processed, or transported. The term facility includes any equipment that may be present along the flow path between a reservoir and its delivery outlets. Facilities may include production wells, injection wells, well tubulars, wellhead equipment, gathering lines, manifolds, pumps, compressors, separators, surface flow lines, steam generation plants, processing plants, and delivery outlets. Examples of facilities include hydrocarbon production fields, polymerization plants, refineries, LNG plants, LNG tanker vessels, and regasification plants, among others.

The term “hydrocarbon” as used herein is an organic compound that primarily includes the elements hydrogen and carbon, although nitrogen, sulfur, oxygen, metals, or any number of other elements may be present in small amounts. As used herein, hydrocarbons generally refer to components found in natural gas, oil, or chemical processing facilities, such as refineries or chemical plants.

The term “natural gas” as used herein refers to a multi-component gas obtained from a crude oil well (associated gas) and/or from a subterranean gas-bearing formation (non-associated gas). The composition and pressure of natural gas can vary significantly. A typical natural gas stream contains methane (CH₄) as a major component, i.e., greater than 50 mol % of the natural gas stream is methane. The natural gas stream can also contain ethane (C₂H₆), higher molecular weight hydrocarbons (e.g., C₃-C₂₀ hydrocarbons), one or more acid gases (e.g., hydrogen sulfide), or any combinations thereof The natural gas can also contain minor amounts of contaminants such as water, nitrogen, iron sulfide, wax, crude oil, or any combinations thereof.

The present disclosure provides methods for differential image quality enhancement for a detection system which includes multiple EM radiation detectors. In one or more embodiments, the detection system may be used to detect one or more chemical species released into the environment. It has been found that, when using detection methods as described in WO2012/134796, there can be variations between the two detectors. The differential image quality enhancement described in the present disclosure increases the contrast in each image while reducing the variation between the images obtained from multiple detectors, thereby improving image quality and helping reduce false positives from background and noise interference. By adjusting the gain, offset, or both of the pixels of each of the plurality of EM radiation detectors based on the intensity values of one or more images from the EM radiation detector configured to detect the emission or absorption of a chemical species (chemical band EM radiation detector) improves the contrast in one or more areas of interest in the images and, therefore, improves the image quality for detection purposes. Embodiments of the present disclosure provide the ability of the detection system to be in motion while detecting a potential leak and provide the ability to detect a potential leak from images collected simultaneously through the use of multiple detectors, allowing for less complex algorithms to be used to analyze for a potential leak. Use of less complex algorithms allows the code to optionally be installed in firmware within components of the collection system, reducing the analysis required by a separate computer system.

Any suitable detection system which includes multiple EM radiation detectors may be used with the differential image quality enhancement method described herein. The multiple EM radiation detectors may include a first EM radiation detector configured to detect a chemical species (e.g., a chemical band EM radiation detector) and a second EM radiation detector configured to provide a reference background (e.g., a reference band EM radiation detector). FIG. 1 illustrates a detection system for detecting a chemical species according to one or more embodiments of the present disclosure. Although one or more embodiments may describe the use of the detection system to detect gaseous hydrocarbon releases such as leaks, the present description is not so limited. The systems and methods of the present disclosure may be used for detecting a chemical species in any application.

Referring to FIG. 1, detection system 10 includes a collection system 55. The collection system 55 including a lens 11 which receives IR rays 21 from a scene to be detected having passed through a gaseous vapor release of a chemical species 25 originating from container 22. The collection system 55 also includes a beam splitter 12 which receives the infrared rays 21 from the lens 11. The beam splitter 12 transmits a first EM radiation beam 18 to a first bandpass filter 13 and reflects a second EM radiation beam 19 to a second bandpass filter 14 reflected off mirror 17. Mirror 17 is optional; however, in the collection system depicted in FIG. 1, mirror 17 is provided so that first EM radiation detector 15 and second EM radiation detector 16 may be arranged parallel to one another (side-by-side). In embodiments without the use of a mirror, the second EM detector would be arranged ninety degrees (90°) from the first EM detector. Optionally, more than one mirror may be used to direct an EM radiation beam within the detection system. The filtered first EM radiation beam 18′ is received by the chemical band EM radiation detector (the first EM radiation detector) 15 and the filtered second EM radiation beam 19′ is received by the reference band EM radiation detector (the second EM radiation detector) 16. The detection system 10 may also include an analysis system 30. The analysis system 30 may include a computer system 35. The computer system 35 includes a display 40, a processor 45, and a non-transitory, computer-readable medium 50. The computer system 35 may be wired to one or more components of the collection system 55 or wirelessly connected to one or more components of the collection system 55. Although shown in FIG. 1 as separate analysis system 30 and computer system 35, in other embodiments, firmware may be installed within a processor; a non-transitory, computer-readable medium; or combinations thereof, within the collection system. In this certain embodiment, the firmware may include the code for directing the processor to perform the differential image quality enhancement as well as any analysis to detect a potential chemical species and other functions. As discussed herein, the use of less complex algorithms allows the code to optionally be installed in firmware within components of the collection system, reducing the analysis required by a separate computer system.

The lens of the detection system may be any suitable lens. For example, the lens may be a fixed focal length lens or a multiple focal length lens (a zoom lens). The lens may be constructed of a material with high transmittance in a desired wavelength range. In one or more embodiments, the lens may be constructed of a silicon material or a germanium material for high transmittance in a mid-IR wavelength range.

The beam splitter may be any suitable beam splitter capable of receiving an EM radiation beam and splitting the EM radiation beam into a plurality of beams, for example at least two or more beams (first and second EM radiation beams or more), or at least three or more beams (first, second, and third EM radiation beams or more), or at least four or more EM radiation beams (first, second, third, and fourth EM radiation beams or more), or at least five or more EM radiation beams (first, second, third, fourth, and fifth EM radiation beams or more). In one or more embodiments, a beam splitter may be used to form two EM radiation beams which may be transmitted to an associated detector in the IR spectrum, and a separate lens may be used to transmit another EM radiation beam to an associated detector in the visible spectrum to infuse any detected chemical species into the visible image for easy human viewing.

In one or more embodiments, the beam splitter may be a broadband splitter or a prism.

In one or more embodiments, the beam splitter may be a broadband beam splitter. The broadband splitter may be constructed from a commercially available material which has approximately a 50% transmittance rate and a 50% reflectance rate. In one or more other embodiments, the broadband splitter may have a higher transmittance rate to one of the detectors. For example, the broadband splitter may have a transmittance rate of 60% or 70% or 80% or more to one of the detectors, while the remaining detectors receive 40% or 30% or 20% or less reflected from the broadband splitter. A broadband splitter transmits and reflects all wavelengths contained in the EM radiation beam received by the broadband beam splitter.

In one or more embodiments, the beam splitter may be a dichroic prism. A dichroic prism transmits the portion of the EM radiation beam with wavelengths below a cutoff wavelength value and reflects the portion of the EM radiation beam with wavelengths above the cutoff wavelength value. The dichroic prism may have a high transmittance rate, for example at least 80% or at least 90% or at least 95% or substantially 100% of the wavelengths below the cutoff wavelength value. The dichroic prism may have a high reflectance rate, for example at least 80% or at least 90% or at least 95% or substantially 100% of the wavelengths above the cutoff wavelength value.

In one or more embodiments, the cutoff value of the dichroic prism may be between a lower limit of the transmittance window of a chemical bandpass filter (e.g., the first bandpass filter) and the upper limit of the transmittance window of a reference bandpass filter (e.g., the second bandpass filter) when the transmittance window of the second bandpass filter is offset (at lower wavelengths) from the transmittance window of the first bandpass filter. In one or more other embodiments, the cutoff wavelength value of the dichroic prism may be between the upper limit of the transmittance window of a chemical bandpass filter (e.g., the first bandpass filter) and the lower limit of the transmittance window of a reference bandpass filter (e.g., the second bandpass filter) when the transmittance window of the second bandpass filter is offset (at greater wavelengths) from the transmittance window of the first bandpass window. For example, the cutoff wavelength value of the dichroic prism may be a value between the upper limit of the first transmittance window of the first bandpass filter and the lower limit of the second transmittance window of the second bandpass filter. As an exemplary embodiment, if the first transmittance window has a wavelength range between 3.25 microns and 3.45 microns (a 200 nm width) and the second transmittance window has a wavelength range between 3.7 microns and 4.1 microns (a 400 nm width), the cutoff wavelength value may be greater than 3.45 microns and less than 3.7 microns, such as 3.5 microns, 3.55 microns, 3.6 microns, or 3.65 microns.

In one or more embodiments, the cutoff wavelength value may be substantially halfway between the offset limits of the first and second transmittance windows. For example, the cutoff wavelength value may be substantially halfway between the upper limit of the first transmittance window and the lower limit of the second transmittance window, or vice versa. In such an arrangement, the dichroic prism transmits the portion of the EM radiation beam having wavelengths less than the cutoff wavelength value and reflects the portion of the EM radiation beam having wavelengths greater than the cutoff wavelength value. In one or more other embodiments, the cutoff wavelength value of the dichroic prism may be greater than the upper limit of the transmittance window of the first bandpass filter or less than the lower limit of the first bandpass filter when the transmittance window for the second bandpass filter overlaps the transmittance window of the first bandpass filter. The cutoff wavelength value of the dichroic prism may be a value as further discussed herein.

In one or more embodiments, a chemical bandpass filter (e.g., a first bandpass filter) may be provided and may have a first transmittance window having a first width. The transmittance window of the first bandpass filter transmits EM radiation within a first EM radiation wavelength range and substantially rejects all other wavelengths outside the first transmittance window. The first EM radiation wavelength range of the first transmittance window may correspond to at least a portion of a wavelength range of absorption or emission by at least one chemical species to be detected. The transmittance window of the first bandpass filter (e.g., a first transmittance window) has a width (e.g., a first width) within an EM radiation wavelength range (e.g., a first EM radiation wavelength range). The transmittance window has a lower limit wavelength value and an upper limit wavelength value for the wavelength range. In one or more embodiments, the transmittance window for the first bandpass filter may have a width of at most 300 nm, at most 250 nm, at most 200 nm, at most 225 nm, at most 150 nm, or at most 100 nm or less. The width of the transmittance window for the first bandpass filter may be in the range of from 25 nm to 300 nm, from 50 nm to 250 nm, or from 75 nm to 200 nm.

The wavelength range of the transmittance window for the first bandpass filter may correspond to at least a portion of the wavelength range of absorption or emission by the at least one chemical species to be detected. In one or more embodiments, the transmittance window for the first bandpass filter is within the EM radiation spectrum, for example the IR spectrum such as the mid-IR spectrum. In one or more embodiments, the transmittance window for the first bandpass filter may have a lower limit of at least 2.5 microns, for example at least 3 microns, at least 3.25 microns, or at least 3.3 microns. In one or more embodiments, the transmittance window for the first bandpass filter may have an upper limit of at most 3.7 microns, for example at most 3.5 microns, at most 3.45 microns, or at most 3.4 microns, for example in the range of from 3 microns to 3.5 microns, from 3.25 microns to 3.45 microns, or from 3.3 microns to 3.4 microns. In one or more embodiments, the chemical species may be a hydrocarbon emitting or absorbing wavelengths within a range of from 3.2 microns to 3.5 microns.

The transmittance window of a chemical bandpass filter (e.g., the first bandpass filter) may transmit any suitable percentage of the first EM radiation beam within the wavelength range of the first transmittance window. In one or more embodiments, the transmittance window of the first bandpass filter may transmit at least 50% of the wavelengths of the first EM radiation beam received by the first bandpass filter within the wavelength range of the transmittance window. The transmittance window may transmit at least 75% or at least 80% or at least 90% or at least 95% or substantially 100% of the first EM radiation beam within the wavelength range of the transmittance window.

In one or more embodiments, a reference bandpass filter (e.g., a second bandpass filter) may be provided and may have a second transmittance window having a second width within an EM radiation wavelength range (e.g., a second EM radiation wavelength range). The second transmittance window passes or transmits EM radiation within a second EM radiation wavelength range and substantially rejects all other wavelengths outside the second transmittance window. The second transmittance window may be of any suitable width. In one or more embodiments, the second transmittance window may have an upper limit wavelength value that is greater than an upper limit wavelength value of the first transmittance window. In one or more embodiments, the width of the second transmittance window may be at least 50 nm, for example at least 100 nm, at least 200 nm, at least 275 nm, at least 300 nm, at least 350 nm, at least 400 nm, or at least 500 nm or more. For example, the width of the transmittance window may be in the range of from 50 nm to 2000 nm, or from 275 nm to 1000 nm, or from 300 nm to 500 nm. In one or more embodiments, the transmittance window for the first bandpass filter may have a width less than, greater than, or equal to the width of the second transmittance window for the second bandpass filter. For example, in certain embodiments, the width of the second transmittance window may be less than the width of the first transmittance window when the second transmittance window has an upper limit wavelength value that is greater than the upper limit wavelength value of the first transmittance window. This decrease in width may be used to reduce the energy level of the filtered EM radiation beam having greater wavelengths (e.g., the filtered second EM radiation beam or the filtered reference band EM radiation beam), also helping to match the energy levels of the filtered first EM radiation beam and the filtered second EM radiation beam.

In one or more embodiments, the transmittance window of a reference bandpass filter (e.g., the second bandpass filter) may partially overlap at the upper limit or the lower limit of the first transmittance window or may completely overlap the transmittance window of a chemical bandpass filter (e.g., the first bandpass filter). For example, the lower limit of the second transmittance window (having an upper limit greater than the upper limit of the first transmittance window) may be less than the upper limit or less than the lower limit of the first transmittance window. In another example, the upper limit of the second transmittance window (having a lower limit less than the lower limit of the first transmittance window) may be greater than the lower limit or greater than the upper limit of the first transmittance window. In one or more embodiments, the second transmittance window may have a lower limit of at least 2 microns, at least 3 microns, at least 3.3 microns or more and an upper limit of at most 5 microns, at most 4.5 microns, or at most 4 microns, for example in the range of from 3 microns to 5 microns, from 3 microns to 4 microns, or from 3.3 microns to 4 microns.

In one or more other embodiments, the transmittance window of a reference bandpass filter (e.g., the second bandpass filter) may be offset from the transmittance window of a chemical bandpass filter (e.g., the first bandpass filter). The amount of offset may vary; however, the offset may be sufficient to enable the EM radiation detectors to distinguish between the background of a scene and the chemical species to be detected within the scene. In other words, the transmittance window of the second EM radiation detector has an offset such that images from the second EM radiation detector include background EM radiation and are devoid of any emission or absorbance from chemical species to be detected. In one or more embodiments, the second transmittance window may be offset from the first transmittance window by at least 100 nm or at least 150 nm or at least 200 nm or more. For example, the lower limit wavelength value of the second transmittance window is greater than the upper limit wavelength value of the first transmittance window, e.g., the lower limit wavelength value of the second transmittance window may be at least 100 nm or at least 150 nm or at least 200 nm or more greater than the upper limit wavelength value of the first transmittance window. For example, the transmittance window for the first bandpass filter may be within the IR spectrum, for example the mid-IR spectrum as discussed herein, and the transmittance window for the second bandpass filter may have a lower limit of at least 3.5 microns or at least 3.7 microns and an upper limit of at most 4.5 microns or at most 4.25 microns or at most 4.1 microns or at most 4 microns, for example in the range of from 3.5 microns to 4.5 microns or from 3.7 microns to 4.1 microns.

The widths and wavelength ranges for the first transmittance window and the second transmittance window (e.g., partial overlap, complete overlap, or offset) of a chemical bandpass filter and a reference bandpass filter may vary by application, equipment and chemical species to be detected, but may at least be sufficient to enable the detectors to distinguish between the background and the chemical species.

The transmittance window of a reference bandpass filter (e.g., the second bandpass filter) may transmit any suitable percentage of the EM radiation beam (e.g., the second EM radiation beam or the reference band EM radiation beam) within the wavelength range of the transmittance window (e.g., the second transmittance window). In one or more embodiments, the transmittance window of a reference bandpass filter (e.g., the second bandpass filter) may transmit a lesser percentage of the EM radiation beam (e.g., the second EM radiation beam) received by the bandpass filter than a chemical bandpass filter (e.g., the first bandpass filter) transmits of the EM radiation beam (e.g., the first EM radiation beam or the chemical band EM radiation beam) through its transmittance window, as described in U.S. Provisional Application Ser. No. 62/015,924, filed Jun. 23, 2014, and titled, “System for Detecting a Chemical Species and Use Thereof”, which is incorporated herein by reference in its entirety. The second transmittance window may transmit less than 50% of the wavelengths within the wavelength range of the transmittance window, for example the second transmittance window may transmit at most 45%, at most 40%, at most 35%, at most 30%, at most 25%, or at most 20% of the second EM radiation beam within the wavelength range of the transmittance window. By decreasing the percentage of the second EM radiation beam transmitted to the detector by lowering the transmittance of the second bandpass filter, the energy level of the filtered second EM radiation beam more closely matches the energy level of the filtered first EM radiation beam, allowing substantially similar integration times and improving image quality between the chemical band and reference band EM radiation detectors.

The first EM radiation detector may receive at least a portion of the first filtered EM radiation beam from the first bandpass filter to detect the filtered first EM radiation beam. The second EM radiation detector may receive at least a portion of the filtered second EM radiation beam from the second bandpass filter to detect the filtered second EM radiation beam. In one or more other embodiments, the first EM radiation beam and the second EM radiation beam may be received by the associated EM radiation detector from the beam splitter without an intervening bandpass filter. In one or more embodiments, the EM radiation beam may be received by the associated EM radiation detector via one or more mirrors which may be used to direct the beam. The EM radiation detectors (e.g., the first and second EM radiation detectors) may be conventional EM radiation detectors configured to detect EM radiation. The detectors include one or more sensors configured to detect EM radiation in the appropriate wavelength range. The one or more sensors are configured to include a plurality of pixels to detect the intensity of the EM radiation received by the detector. For example, the sensors may include a focal plane array, a charge-coupled device (CCD), a complementary metal-oxide-semiconductor (CMOS), and any combinations thereof In one or more embodiments, the one or more sensors of the EM radiation detectors may be included within a camera. In one or more embodiments, the camera may include the components described herein within the collection system. The collection system or camera may also include a readout integrated circuit (ROIC), video interface board, a field programmable gate array (FPGA), and any combinations thereof. The collection system or camera may also include a cooler such as a cryogenic cooler, for example an integrated dewar cooler assembly (IDCA), configured to cool the sensors in the EM radiation detectors. The EM radiation detectors are configured to communicate with the analysis system of the detection system. The communication components may provide wired or wireless communication.

In one or more embodiments, the EM radiation detectors may be configured to detect mid-IR EM radiation. The mid-IR EM radiation detector may include one or more focal plane arrays which may include a mercury-cadmium-telluride focal plane array, an indium-Antimonide focal plane array, an indium-gallium-arsenide focal plane array, a vanadium oxide focal plane array, and any combinations thereof. The mid-IR EM radiation detectors may include at least one cooler configured to cool the detector. Examples of gas-detecting IR cameras that are commercially available include GF320 and GF306 cameras manufactured by FLIR Systems, Inc., and EYE-C-GAS® camera manufactured by Opgal.

In one or more other embodiments, the multiple EM radiation detectors may include an EM radiation detector configured to detect visible or near-IR EM radiation. The visible or near-IR EM radiation detector may be a CCD, a CMOS, and any combinations thereof.

In one or more embodiments, the EM radiation beam forms an image on the associated EM radiation detector. The EM radiation detector may be configured to generate a single image, a plurality of still images, or a video of sequential images. Any number of EM radiation detectors may be used in a detection system, depending on the wavelength ranges to be detected. Thus, the multiple EM radiation detectors may be the same or different.

In one or more embodiments, an EM radiation detector of the detection system may be mounted on a pole with at least two degrees of freedom, such as panning and tilting, under the manual control or the automatic control of the system. Movement of an EM radiation detector may include 360 degree coverage. Several EM radiation detectors may be positioned around the perimeter of a plant to give 100% coverage of the facility. An autonomous system can provide continuous plant surveillance. The system can continue monitoring the facility 24 hours a day, seven days a week, and 365 days per year with minimal downtime. Downtime may mainly be the result of performing routine maintenance on the system and may be compensated for by redundancy, e.g., directing one or more other EM radiation detectors at an area whose EM radiation detectors are being serviced. The overall system cost may be kept low while keeping the false alarm rate low and still being able to detect small or early hydrocarbon leaks, e.g., plumes with about 20% LEL at a distance of 150 meters subject to environmental conditions.

The detection system is not limited to pole-mounted EM radiation detectors. In other embodiments, the EM radiation detectors may be attached to autonomous mobile platforms, placed on conveniently located towers, or suspended from cables or balloons. The detection system may also be integrated into mobile robots, which are either autonomous or steered by an operator.

In one or more embodiments, the detection system may be calibrated. Any suitable calibration method may be used. In one or more embodiments, a differential calibration method may be used, as described in U.S. Provisional Application Ser. No. 62/015,926, filed Jun. 23, 2014, and titled, “Methods for Calibrating Multiple Detector System”, which is incorporated herein by reference in its entirety. In the differential calibration method, a multiple detector system may be calibrated by matching individual pixels of an image from a chemical band EM radiation detector (e.g., the first EM radiation detector) and individual pixels of an image from a reference band EM radiation detector (e.g., the second EM radiation detector) to the average intensity of the pixels of the chemical band EM radiation detector. By matching the individual pixel intensity of the reference band EM radiation detector (e.g., the second EM radiation detector) to the average intensity of the chemical band EM radiation detector (e.g., the first EM radiation detector), an improved image quality can be obtained when using the detection system to determine the presence of a chemical species. This differential calibration method can decrease variability between the multiple detectors, increasing the ability to detect the presence of a chemical species that might otherwise be masked by the intensity response variability between detectors. For example, the differences in intensity between the chemical band EM radiation detector and the reference band EM radiation detector from independent calibration could be larger than the signal produced by the presence of a chemical species. This differential calibration method also improves image quality and helps reduce false positives from background and noise interference

A calibration EM radiation beam may be generated using a calibration EM radiation source. The calibration EM radiation source may be any suitable source capable of providing a substantially uniform source of EM radiation for calibrating an EM radiation detector. The calibration EM radiation source is generated across the entire field of view of the EM radiation detectors. For example, the calibration EM radiation source may be a blackbody. The calibration EM radiation beam is generated by the source at a first temperature. The first temperature may be in the range of from minus 50° Celsius (C) to 250° C. or from 0° C. to 175° C. or from 5° C. to 150° C. or from 20° C. to 100° C.

At least a portion of the calibration EM radiation beam may be detected with a chemical band EM radiation detector (e.g., the first EM radiation detector). The calibration EM radiation beam may be provided to the multiple EM radiation detectors by passing through a beam splitter to provide a plurality of calibration EM radiation beams (e.g., a first calibration EM radiation beam and a second calibration EM radiation beam). At least a portion of the first calibration EM radiation beam may pass through a chemical bandpass filter (e.g., the first bandpass filter) to generate a filtered first calibration EM radiation beam. At least a portion of the second calibration EM radiation beam may pass through a reference bandpass filter (e.g., the second bandpass filter) to generate a filtered second calibration EM radiation beam. An average intensity (mean intensity) value of the calibration EM radiation beam detected by the first EM radiation detector may be obtained. The average intensity value may be determined by calculating the average intensity of the plurality of pixels contained within the first EM radiation detector.

One or more of the plurality of pixels of the image from the first EM radiation detector may be adjusted to decrease the difference between the intensity of an individual pixel and the average intensity value of the first EM radiation detector. One or more of the plurality of pixels of the image from the second EM radiation detector may be adjusted to decrease the difference between the intensity of an individual pixel and the average intensity value of the first EM radiation detector. The intensity of an individual pixel may be adjusted by modifying the gain, offset, and any combinations thereof.

In one or more embodiments, after the calibration adjustment of the pixels of the first EM radiation detector, the difference between an individual pixel intensity and the average pixel intensity of the first EM radiation detector may be within at most 5% of the dynamic range of the sensor of the first EM radiation detector, for example within 1% or within 0.5% or within 0.1%.

In one or more embodiments, after the calibration adjustment of the pixels of the second EM radiation detector, the difference between an individual pixel intensity and the average pixel intensity of the first EM radiation detector may be within at most 5% of the dynamic range of the sensor of the second EM radiation detector, for example within 1% or within 0.5% or within 0.1%. In one or more embodiments, the intensity of a majority (more than 50%) of the plurality of pixels of the first EM radiation detector and a majority (more than 50%) of the plurality of pixels of the second EM radiation detector are substantially the same as the average intensity value for the first EM radiation detector after calibration adjustment.

If additional calibration is desired for additional temperatures, the method reverts back to the beginning and generates a calibration EM radiation beam at one or more additional temperatures. This may be repeated for any number of desired calibration temperatures. The additional calibration temperatures may be within the ranges discussed herein for the first calibration temperature but have different values. If no additional calibration temperatures are desired, the differential calibration method may end. Additionally, if more than two detectors are used in the detection system, the method may also include additional steps to calibrate the additional detectors, for example additional detectors similar to the first EM radiation detector may be calibrated by adjusting one or more pixels of the additional EM radiation detector to decrease the difference between the intensity of an individual pixel and the average intensity value of the first EM radiation detector, and additional detectors which are not similar to the first EM radiation detector, such as a visible EM radiation detector, may be calibrated using any appropriate calibration method for such detector.

In one or more embodiments, the code of the analysis system may additionally be configured to direct the processor to: obtain an average intensity value of a plurality of pixels from a chemical band EM radiation detector (e.g., the first EM radiation detector) detecting the calibration EM radiation beam; adjust one or more pixels of an image from the first EM radiation detector to decrease the difference between the intensity of an individual pixel and the average intensity value of the first EM radiation detector; and adjust one or more pixels of a reference band EM radiation detector (e.g., the second EM radiation detector) detecting the calibration EM radiation beam to decrease the difference between the intensity of an individual pixel and the average intensity value of the first EM radiation detector.

According to the present disclosure, the detection system includes an analysis system configured to analyze the output from the multiple EM radiation detectors. The analysis system is configured to provide differential image quality enhancement for the multiple detectors. The detection system may include an analysis system configured to analyze the output from the EM radiation detectors. The analysis system may include a processor and non-transitory, computer-readable medium. The processor, non-transitory, computer-readable medium or combinations thereof may comprise code. The analysis system may also include a display. The analysis system may also include a graphical processing unit (GPU). The code is configured to direct the processor to: identify an image from one or more images from a chemical band EM radiation detector (e.g., the first EM radiation detector), the image including a plurality of pixels, each pixel having an associated intensity value; identify an image from one or more images from a reference band EM radiation detector (e.g., the second EM radiation detector), the image including a plurality of pixels, each pixel having an associated intensity value; adjust one or more intensity values of the plurality of pixels of the image from the second EM radiation detector based on one or more intensity value parameters of the image from the first EM radiation detector; and adjust one or more intensity values of the plurality of pixels of the image from the first EM radiation detector based on one or more intensity value parameters of the image from the first EM radiation detector. The differential image quality enhancement may be applied to the images automatically or through manual operator input. The adjustment values may be continually updated or periodically updated after a set number of images have been generated, for example no more than 50 successive images have been generated, or no more than 25 successive images have been generated, or no more than 15 successive images have been generated. With periodic updating, the determined adjustment values for the pixels may be applied to the pixels of subsequent images to enhance the image quality of those images until the adjustment values are updated. The code of the analysis system may be further configured to direct the processor to perform various other functions of methods described herein.

The method for differential image quality enhancement includes obtaining an image from one or more images from a chemical band EM radiation detector (e.g., the first EM radiation detector). The chemical band image includes a plurality of pixels, each pixel having an associated intensity value. An image from one or more images from a reference band EM radiation detector (e.g., the second EM radiation detector) is obtained. The reference band image includes a plurality of pixels, each pixel having an associated intensity value. One or more intensity values of the plurality of pixels of the image from the second EM radiation detector are adjusted based on one or more intensity value parameters of the image from the first EM radiation detector. The adjusting of the one or more intensity values of the plurality of pixels may include adjusting the gain, offset, and combinations thereof based on one or more intensity value parameters of the image from the first EM radiation detector. The intensity value parameters may include a maximum intensity value for the plurality of pixels (I_(max)), a minimum intensity value for the plurality of pixels (I_(min)), and combinations thereof.

In one or more embodiments, the maximum intensity value, the minimum intensity value, and combinations thereof may be determined based on a subset of the plurality of pixels from an image from a chemical band EM radiation detector corresponding to an area of interest. A plurality of areas of interest may be identified in images relevant for applying the differential image quality enhancement. A maximum intensity value, a minimum intensity value, and combinations thereof may be determined for each of the areas of interest based on corresponding subsets of the plurality of pixels from an image from a chemical band EM radiation detector. Corresponding subsets of the plurality of pixels from an image from a reference band EM radiation detector may be adjusted based on the maximum intensity values, minimum intensity values, and combinations thereof for each of the areas of interest of the image from the chemical band EM radiation detector.

An EM detector has an input frame having a given range of input frame intensities and an output frame having a given range of output frame intensities. The range of intensities for the output frame may be less than the range for the input frame, requiring an adjustment of the image. Adjustment of the gain and/or offset of the pixels of an image, whether from a chemical band EM radiation detector or a reference band EM radiation detector, may be based on the intensity values of the pixels of the chemical band image (e.g., the first plurality of pixels). The plurality of pixels in the images from the chemical band EM radiation detector and the plurality of pixels in the images from the reference band EM radiation detector may substantially correspond temporally, spatially, or both. The gain and offset may be adjusted for individual pixels from the multiple detectors based on the intensity values of the plurality of pixels of the chemical band image. Pixels from subsequent images from the chemical band EM radiation detector and the reference band EM radiation detector (e.g., the first EM radiation detector and the second EM radiation detector) may also be adjusted based on the values for gain and/or offset determined from the differential image quality enhancement analysis of one or more prior images. In one or more embodiments, the differential image quality enhancement method may be applied periodically to update values for gain and/or offset to be applied to the pixels of the images. The frequency of which the method may be applied may be after a set number of images have been generated, for example the periodic adjustment occurs after no more than 50 successive images have been generated, or no more than 25 successive images have been generated, or no more than 15 successive images have been generated.

The adjustment of the gain of a pixel depends on the dynamic range of the input frame and the dynamic range of the output frame. The maximum intensity measured by the input frame and the minimum intensity measured by the input frame are determined, and then a differential dynamic range representing the difference between the maximum intensity and the minimum intensity is determined. The intensities of the pixels in the output frame are adjusted based on the differential dynamic range of intensities determined from the input frame. For example, the input frame of the first and second EM radiation detectors may be a video frame having a 14-bit depth (2¹⁴ for a 0 to 16384 dynamic intensity range) and an output frame of an 8-bit depth (2⁸ for a 0 to 256 dynamic intensity range). If the first plurality of pixels of the image from the first EM radiation detector has a maximum intensity of 7100 and a minimum intensity of 4000 in the input frame, the differential dynamic range is 3100. The output frame has a dynamic intensity range of 0 to 256, so the differential dynamic range of 3100 may be scaled to the 256 dynamic range of the output frame (gain control). This technique may also be applied separately for a plurality of areas of interest within the chemical band image using the intensity values within a subset of pixels for each area of interest. The differential gain control may be achieved using any suitable method, for example the gain to be applied to the pixel intensities in the output frame for a 256 dynamic intensity range may be determined by the following:

Gain=256/(I_(max)−I_(min)), where I_(max) is the maximum intensity pixel from the first plurality of pixels of the first EM radiation detector and I_(min) is the minimum intensity pixel from the first plurality of pixels of the first EM radiation detector.

The offset may be adjusted by adding or subtracting a quantity to the intensity of a pixel. Offset may be determined using any suitable method, for example the offset may be determined with respect to the example discussed above by the following:

Offset=I _(min)*(256/(I _(max) −I _(min))).

The intensity of the pixel in the output frame may be determined by the following: I_(out)=(I*Gain)−Offset, where I is the intensity of the pixel in the input frame and I_(out) is the intensity of the pixel in the output frame. Therefore, in the example, for a pixel of the second plurality of pixels from the second EM radiation detector which has an intensity in the output frame of 6435, the gain would be 0.08 and the offset would be 330 so that the intensity in the output frame (I_(out)) would be 185. In this example, the other pixels of the second plurality and first plurality are likewise adjusted to achieve corresponding output frame intensity levels.

In one or more other embodiments, a histogram of the intensity values of the pixels may be generated. FIGS. 2A and 2B illustrate histograms of a first plurality of pixels 205 of an image from a chemical band EM radiation detector (e.g., a first EM radiation detector) configured to detect wavelengths in the range for the emission or absorption of a chemical species and a second plurality of pixels 210 of an image from a reference band EM radiation detector (e.g., a second EM radiation detector) configured to detect wavelengths outside of the range of wavelengths of the emission or absorption of a chemical species. As shown in FIG. 2A, substantially all of the first plurality of pixels have intensities within two peaks 215, 220. The intensity is represented along the x-axis and the number of pixels having the same intensity value is represented along the y-axis. The second plurality of pixels from the second EM radiation detector similarly has substantially all of the pixels having intensities within two peaks 225, 230; however, peaks 225, 230 are shifted to the right of peaks 215, 220. An offset may be applied to either the first plurality of pixels or the second plurality of pixels to adjust the intensity values of the plurality of pixels such that the intensity value distribution of the first plurality of pixels substantially aligns with the intensity value distribution of the second plurality of pixels (histogram matching method). As shown in FIG. 2B, peaks 215, 225 and 220, 230 are substantially aligned after applying the offset. The intensity values of the second plurality of pixels and first plurality of pixels may also be adjusted by adjusting the gain. The gain may be determined as discussed herein.

In one or more embodiments, one or more subsets of the first plurality of pixels having similar intensities may be identified as one or more areas of interest (e.g., a first area of interest, a second area of interest, etc.) within the histogram. The corresponding pixels of the second plurality of pixels are also identified. The one or more areas of interest may each consist of one or more pixels having intensity values within a defined range of intensity values. For example, referring to FIGS. 2A and 2B, a first area of interest may include both peaks 215, 220 for the first plurality of pixels and peaks 225, 230 for the second plurality of pixels or a first area of interest may include peak 215 for the first plurality of pixels and peak 225 for the second plurality of pixels, and a second area of interest may include peak 220 for the first plurality of pixels and peak 230 for the second plurality of pixels.

In one or more other embodiments, the images from the first EM radiation detector and the second EM radiation detector may be separated into a plurality of areas of interest. The plurality of areas of interest may include a majority (more than 50%) of the area of the images.

In one or more other embodiments, the pixels within one or more areas of interest (e.g., a first area of interest, a second area of interest, etc.) may be determined by identifying at least one of the first plurality of pixels along with one or more adjacent pixels to form a subset of pixels defining an area of interest. The corresponding pixels of the second plurality of pixels may be identified for each subset and a histogram of the intensity values of the pixels within each of the areas of interest may be generated. The pixels of the first plurality of pixels and the second plurality of pixels within the one or more areas of interest may be substantially aligned temporally, spatially, or both. An offset may be applied based on the intensity values within each of the one or more areas of interest to substantially align the intensity distributions. The gain may also be determined based on the intensity values within each of the one or more areas of interest. The one or more areas of interest may be of any suitable size, for example at least a 10×10 grid of adjacent pixels, at least a 20×20 grid, etc.

Identifying the pixels for a subset may be determined by the differences between the intensity values of the first plurality of pixels and the corresponding second plurality of pixels and applying an image enhancement threshold value to the determined differences in intensity values. One or more pixels satisfying the image enhancement threshold value may be included in a subset for one of the one or more areas of interest. The image enhancement threshold value may be similar to the resultant threshold values described herein.

In one or more embodiments, the differences between intensity values between the plurality of pixels (e.g., n₁ to n_(z)) of an image from a chemical band EM radiation detector (e.g., the first EM radiation detector) and the corresponding plurality of pixels (e.g., m₁ to m_(z)) of an image from a reference band EM radiation detector (e.g., the second EM radiation detector) may be determined by calculating the ratio of the intensity values between the corresponding pixels of the images. For example, the differences between intensity values between the first plurality of pixels (e.g., n₁ to n_(z)) of the first image and the corresponding second plurality of pixels from the second image may be determined by calculating the ratio of the intensity values between each of the first plurality of pixels of the first image and the corresponding second plurality of pixels of the second image.

In one or more other embodiments, the differences between intensity values between the plurality of pixels (e.g., n₁ to n_(z)) of an image from a chemical band EM radiation detector (e.g., the first EM radiation detector) and the corresponding plurality of pixels (e.g., m₁ to m_(z)) of an image from a reference band EM radiation detector (e.g., the second EM radiation detector) may be determined by subtracting the intensity values between the corresponding pixels of the images. For example, the differences between intensity values of the first plurality of pixels and the corresponding second plurality of pixels is determined by subtracting the intensity values between each of the first plurality of pixels (e.g., n₁ to n_(z)) of the first image and the corresponding second plurality of pixels of the second image. The difference in intensity values may be represented by both the sign and the magnitude or may be absolute values.

An offset may be applied to either the first plurality of pixels or the second plurality of pixels within each of the one or more areas of interest to adjust the intensity values of the plurality of pixels such that the intensity value distribution of the first plurality of pixels substantially aligns with the intensity value distribution of the second plurality of pixels within the area of interest. The intensity values of the second plurality of pixels and first plurality of pixels may also be adjusted by adjusting the gain based on the intensity value parameters of the pixels of the chemical band image. The gain may be determined based on the intensity value parameters of the first plurality of pixels as a whole or the intensity value parameters of the subset of the first plurality of pixels within each of the areas of interest. The subsets of pixels of the first plurality of pixels and the second plurality of pixels within the one or more areas of interest may be substantially aligned temporally, spatially, or both.

In one or more other embodiments, differential non-uniformity correction (DNUC) methods may be used in the systems and methods for differential image quality enhancement. As discussed further herein, the differential non-uniformity correction methods use pixel values from images from multiple detectors to identify one or more areas of interest to determine the gain and offset adjustments for the images as well as subsequent images. Such method may include identifying one or more areas of interest for differential image quality enhancement by providing a calibration EM radiation source and obtaining a baseline image (e.g., a first baseline image) from a chemical band EM radiation detector (e.g., the first EM radiation detector), the detector detecting at least a portion of the calibration EM radiation source at a first temperature. The first baseline image includes a plurality of pixels, each pixel having an associated intensity value. A baseline image (e.g., a second baseline image) may be obtained from a reference band EM radiation detector (e.g., the second EM radiation detector), the detector detecting at least a portion of the calibration EM radiation source at a first temperature. The second baseline image includes a plurality of pixels, each pixel having an associated intensity value. Baseline information of the difference in intensity values between each of the plurality of pixels of the first baseline image and the plurality of pixels of the second baseline image may be determined. The plurality of pixels of the first baseline image and the plurality of pixels of the second baseline image may be substantially aligned temporally, spatially and combinations thereof. A table of the baseline information may be generated (e.g., a first table). The table includes the differences in intensity values between the plurality of pixels of the first baseline image and the plurality of pixels of the second baseline image. Upon further generating images of a scene in which a chemical species may be detected, the difference in intensity values between the first plurality of pixels of the image of the scene (e.g., the first image) from a chemical band EM radiation detector (e.g., the first EM radiation detector) and the second plurality of pixels of the image of the scene (e.g., the second image) from a reference band EM radiation detector (e.g., the second EM radiation detector) may be determined (e.g., a resultant image). The first plurality of pixels and the second plurality of pixels may be substantially aligned temporally, spatially, and combinations thereof. The baseline information of the table may then be applied to the difference in intensity values between the first plurality of pixels and the corresponding second plurality of pixels of the images of the scene to be detected. Pixels proximate one another representing non-zero values after application of the table values may be grouped to form the one or more areas of interest. The intensity values of the first plurality of pixels and the second plurality of pixels may also be adjusted according to the embodiments discussed herein with respect to gain and/or offset within each of the one or more areas of interest based on one or more intensity value parameters of the first image of the scene.

In one or more other embodiments, the DNUC method may be based on a table from images of a scene for detection instead of images from a calibration EM radiation source. In this embodiment, an image (e.g., a first baseline image) may be obtained from a chemical band EM radiation detector (e.g., the first EM radiation detector) of the scene for detection (in the absence of a chemical species to be detected), and an image (e.g., a second baseline image) may be obtained from a reference band EM radiation detector (e.g., the second EM radiation detector) of the scene for detection. Baseline information of the difference in intensity values between each of the first plurality of pixels of the first baseline image from the first EM radiation detector and the second plurality of pixels of the second baseline image from the second EM radiation detector may be obtained. The first plurality of pixels and the second plurality of pixels may be substantially aligned temporally, spatially, and combinations thereof A table (e.g., a first table) of the baseline information may be generated. The table includes the differences in intensity values between the first plurality of pixels and the second plurality of pixels. Subsequent additional images from the first and second EM radiation detectors may be generated. The additional images from the first EM radiation detector (e.g., a first image, a third image, a fifth image, etc.) and additional images from the second EM radiation detector (e.g., a second image, a fourth image, a sixth image, etc.) may be obtained. The difference in intensity values between corresponding images from the first EM radiation detector and the second EM radiation detector may be determined, and the baseline information of the table may then be applied to the difference in intensity values of the plurality of pixels. Pixels proximate one another representing non-zero values after application of the table values may be identified and may be grouped to form the one or more areas of interest. The intensity values of the plurality of pixels from images from the chemical band EM radiation detector and images from the reference band EM radiation detector may be adjusted according to the embodiments discussed herein with respect to gain and/or offset within each of the one or more areas of interest based on one or more intensity value parameters of the image from the chemical band EM radiation detector.

For example, a first baseline image from the first EM radiation detector and a second baseline image from the second EM radiation detector of a scene for detection may be used to generate a first table of baseline information. A first image from the first EM radiation detector of a scene for detection of a chemical species and a second image from the second EM radiation detector may be obtained. The first image includes a first plurality of pixels and the second image includes a second plurality of pixels. Each of the first plurality of pixels and the second plurality of pixels are temporally and spatially aligned. The baseline information in the first table is applied to the difference in intensity values between each of the temporally and spatially aligned first plurality of pixels and second plurality of pixels. The pixels proximate one another representing non-zero values are identified and may be grouped into a subset to form one of the one or more areas of interest. Within each of the identified one or more areas of interest, the intensity values of the subset of the first plurality of pixels of the first detection image and the second plurality of pixels of the second detection image are adjusted based on one or more intensity value parameters of the first detection image for each of the one or more areas of interest, according to embodiments discussed herein. A third detection image from the first EM radiation detector and a fourth detection image from the second EM radiation detector may be obtained. The third detection image includes a third plurality of pixels and the fourth detection image includes a fourth plurality of pixels. Each of the third plurality of pixels and the fourth plurality of pixels are temporally and spatially aligned. The baseline information in the first table is applied to the difference in intensity values between each of the temporally and spatially aligned third plurality of pixels and fourth plurality of pixels. The pixels proximate one another representing non-zero values are identified and may be grouped into a subset to form one of the one or more areas of interest. Within each of the identified one or more areas of interest, the intensity values of the third plurality of pixels of the third detection image and the fourth plurality of pixels of the fourth detection image are adjusted based on one or more intensity value parameters of the third detection image for each of the one or more areas of interest, according to embodiments discussed herein.

In one or more other embodiments, differential high sensitivity mode (DHSM) methods may be used in the systems and methods for differential image quality enhancement. As discussed further herein, the DHSM methods use pixel intensity values from images from multiple detectors to identify one or more areas of interest to determine the gain and/or offset adjustments for the images as well as subsequent images. In one or more embodiments, the DHSM method determines a moving average of the differences in intensity values of the corresponding pixels and uses the moving average as the baseline information in the table. The baseline information is applied to the differences in intensity values of subsequent images from the chemical band EM radiation detector (e.g., the first EM radiation detector) and the reference band EM radiation detector (e.g., the second EM radiation detector). Pixels proximate one another representing non-zero values after application of the table values may be identified and may be grouped to form the one or more areas of interest. Within each of the identified one or more areas of interest, the intensity values of the plurality of pixels of a subsequent image from the first EM radiation detector and the second EM radiation detector may also be adjusted based on one or more intensity value parameters of the corresponding subsequent image from the first EM radiation detector, according to the embodiments discussed herein. After the identification of the one or more areas of interest, the table may be updated to include the subsequent images in the moving average of the baseline information.

For example, a first, third, and fifth image may be obtained from the first EM radiation detector and a second, fourth, and sixth image may be obtained from the second EM radiation detector. The first and second images are substantially aligned temporally and spatially. The third and fourth images are substantially aligned temporally and spatially. The fifth and sixth images are substantially aligned temporally and spatially. The first, second, third and fourth images are used as baseline images for determining the moving averages for the baseline information. The moving average is determined based on the differences in intensity values of the first plurality of pixels of the first image and the second plurality of pixels of the second image and the third plurality of pixels of the third image and the fourth plurality of pixels of the fourth image. The fifth and sixth images are used as detection images for determining the presence of a chemical species. The difference in intensity values between the fifth plurality of pixels of the fifth image and the sixth plurality of pixels of the sixth image is determined and the baseline information is applied. Pixels proximate one another representing non-zero values after application of the table values are identified and may be grouped into a subset to form one of the one or more areas of interest. Within each of the identified one or more areas of interest, the intensity values of the plurality of pixels of the subsequent image (e.g., fifth image) from the first EM radiation detector and the subsequent image (e.g., sixth image) from the second EM radiation detector may be adjusted based on one or more intensity value parameters of the corresponding subsequent image (e.g., fifth image) from the first EM radiation detector, according to the embodiments discussed herein. After identification of the one or more areas of interest, the moving averages of the baseline information are updated to include the difference in intensity values of the fifth and sixth images, etc. This may continue for any number of subsequent images from the first and second EM radiation detectors. In one or more embodiments, the moving averages for the baseline information of the differences in intensity values may be determined based on at least 10 sets of previous images from the first and second EM radiation detectors, or at least 25 sets of images, or at least 100 sets of images.

In one or more embodiments of the DNUC and DHSM methods, a baseline threshold value may be applied to the pixels representing non-zero values after applying the baseline information from the table to the difference in intensity values of the pixels of the images. Applying a baseline threshold value may eliminate pixels being identified for the one or more areas of interest due to noise and interference variations.

FIG. 3 illustrates a flow chart for a method of differential image quality enhancement for a detection system including multiple EM radiation detectors in accordance with one or more embodiments of the present disclosure. The differential image quality enhancement method 300 begins at block 305. At block 310, an image is obtained from a chemical band EM radiation detector (e.g., the first EM radiation detector). At block 315, an image is obtained from a reference band EM radiation detector (e.g., the second EM radiation detector). At block 320, one or more intensity values of a plurality of pixels of the image from the reference band EM radiation detector are adjusted (gain and/or offset) based on one or more intensity value parameters from the image from the chemical band EM radiation detector, according to one or more embodiments described herein. At block 325, one or more intensity values of a plurality of pixels of the image from the chemical band EM radiation detector are adjusted (gain and/or offset) based on the one or more intensity value parameters from the chemical band image, according to one or more embodiments described herein. At block 330, the adjustment values applied at block 320, 325 may be applied to subsequent images from the chemical band and reference band EM radiation detectors (e.g., the first and second EM radiation detectors). The adjustment values from block 320, 325 may be periodically updated by returning to block 310 of the differential image quality enhancement process. If no further adjustments are desired, the enhancement process ends at block 335.

The detection system may be used in any facility that has hydrocarbons, or other detectable chemical species, present. Examples of such facilities include LNG plants, oil and gas wellhead operations, offshore platforms, transport pipelines, ships, trucks, refineries, and chemical plants, as described herein. As noted, the chemical plume may be a gaseous hydrocarbon or an oil slick on a surface of water, such as around an offshore platform, tanker, off-loading platform, and the like.

If a positive identification of a leak is made, the detection system may locate the leak and activate an alarm, alerting an operator to send a response team to the site of the leak. The response team may confirm the presence of the leak and effectuate repairs. In one or more embodiments, the hydrocarbon leak may be shown as a false color image for easier operator interpretation. Further, the system may have zoom capability to assist the operator when doing a leak investigation in a manual mode. In one or more embodiments, the EM radiation detectors may be able to be operated in both the automatic and manual modes. Thus, in the event of an alarm, an operator may be able to take control of the EM radiation detectors to do further investigation.

In one or more embodiments, the system may be configured to work over a broad temperature range, including warm and cold temperatures such as a hot, tropical, or desert environment or a cold, arctic environment. Further, the detection system may be adapted to function in the day or night and at a variety of temperatures, for example ranging from minus 50° C. to 100° C. The detection system may also be configured to operate under other environmental interferences such as in fog, rain, or sandstorms. In one or more embodiments, the detection system may detect hydrocarbons such as methane, ethane, propane, butane, pentane, hexane, heptane, octane, ethylene, propylene, isoprene, benzene, ethyl benzene, toluene, xylene, and methyl ethyl ketone, among others. The detection system may be additionally or alternatively configured to detect other chemical species which are capable of being imaged.

The present detection systems may utilize ambient EM radiation energy for the detection, but may also utilize artificial illumination sources for the EM radiation to detect a chemical species. In one or more embodiments, an EM radiation source, e.g., a light source, may be used to illuminate the environment. For example, an IR laser may be used to illuminate an area of interest for leak confirmation. The light source may be useful in conditions in which the contrast between a plume and the background may not be sufficient to distinguish the chemical species. The light source may be powered, activated, or moved using a light source control in communication with the present system.

The detection system is not limited to the detection of chemical plumes, but may also provide other functionality. For example, in one or more embodiments, the detection system may be used to detect an event other than the release of a chemical species such as to monitor specific equipment, such as furnaces, reactors, compressors, and the like, looking for such problems as hot spots, maldistribution, hot motors, and the like. Further, the detection system may provide fence-line monitoring for security purposes and monitoring of fugitive emissions into the environment from equipment.

In one or more embodiments, the detection system may include a master clock. The EM radiation detectors may be configured to receive a signal from the master clock to synchronize the frame rate and/or integration time of the EM radiation detectors. The synchronization may be accomplished through the clocking circuit in the ROIC of the EM radiation detectors. Such synchronization may improve the temporal alignment of the images. In one or more embodiments, one of the EM radiation detectors may be designated as the master and generates a clock signal to the other EM radiation detectors, which are designated as a slave to the master EM radiation detector and receive the clock signal. The resulting frame rate of each EM radiation detector is then synchronized to the shared clock signal.

In one or more embodiments, images from the chemical band EM radiation detector and the reference band EM radiation detector may be spatially aligned. Images may be spatially aligned using any suitable method. In one or more embodiments, spatial alignment may be accomplished by keeping the EM radiation detectors stationary. In one or more other embodiments, the images may be registered such that a software program may provide spatial alignment using a software registration method which uses features in the image for alignment. In one or more embodiments, the EM radiation detectors may be configured to be substantially spatially aligned. In one or more embodiments, the EM radiation detectors may be both substantially temporally and spatially aligned. The images may be registered such that the images may be substantially matched in spatial alignment pixel by pixel, temporal alignment, or the combination thereof

In one or more embodiments, the integration times of the EM radiation detectors may be substantially the same or different. In one or more embodiments, the integration time of a chemical band EM radiation detector (e.g., the first EM radiation detector) may be substantially the same as the integration time of a reference band EM radiation detector (e.g., the second EM radiation detector) and the transmittance of a reference bandpass filter (e.g., the second bandpass filter) may be less than a chemical bandpass filter (e.g., the first bandpass filter) as described herein. By decreasing the transmittance of a reference bandpass filter, the energy level of the resulting filtered EM radiation beam can more closely match the energy level of a filtered chemical band EM radiation beam, assisting in matching the resulting intensity levels between the chemical band EM radiation detector and the reference band EM radiation detector.

In one or more other embodiments, the integration time of a reference band EM radiation detector (e.g., the second EM radiation detector) may be different than the chemical band EM radiation detector (e.g., the first EM radiation detector) to assist in matching the energy level between the reference band EM radiation detector and the chemical band EM radiation detector within a single clock cycle, assisting in matching the resulting intensity levels between the chemical band EM radiation detector and the reference band EM radiation detector. For example, the integration time of a reference band EM radiation detector (e.g., the second EM radiation detector) may be less than the integration time of a chemical band EM radiation detector (e.g., the first EM radiation detector). In one or more embodiments, the integration time of the second EM radiation detector may be at most 95% of the integration time of the first EM radiation detector, for example at most 75%, at most 50%, at most 40%, or at most 33% of the integration time of the first EM radiation detector. When the energy level of the reference band EM radiation beam is greater than the chemical band EM radiation beam, decreasing the integration time of a reference band EM radiation detector (e.g., the second EM radiation detector) can more closely match the energy levels of the EM radiation beams, assisting in matching the resulting intensity levels between the chemical band EM radiation detector and the reference band EM radiation detector. Conversely, when the energy level of the chemical band EM radiation beam is greater than the reference band EM radiation beam, decreasing the integration time of a chemical band EM radiation detector (e.g., the first EM radiation detector) can more closely match the energy levels of the EM radiation beams. However, significant differences in integration times can result in image blurring due to the time differences, reducing image quality.

In one or more embodiments, the code of the analysis system may additionally be configured to direct the processor to: identify at least one image of a scene to be detected by a chemical band EM radiation detector (e.g., the first EM radiation detector) and at least one image of a scene to be detected by a reference band EM radiation detector (e.g., the second EM radiation detector), the images having been subjected to the differential image quality enhancement methods described herein; compare the image from the first EM radiation detector with the image from the second EM radiation detector; and determine if a difference between the images represents the detection of a chemical species.

In one or more embodiments, determining the presence of at least one chemical species may include generating a resultant image. The resultant image may be based, at least in part, on the data generated by the first and second EM radiation detectors, for example a single image, a series of still images, or video images generated by the detectors. The presence or absence of a chemical species may be determined based, at least in part, on any differences existing in the resultant image. Subsequent resultant images may be generated to provide greater accuracy in the identification of a potential release, decreasing false positives resulting from background and noise interference.

In one or more embodiments, a preliminary threshold criteria may be applied to the pixels from images from a chemical band EM radiation detector and images from a reference band EM radiation detector to remove pixels satisfying the preliminary threshold criteria from further consideration. The preliminary threshold criteria may include at least a first preliminary threshold value and a second preliminary threshold value. Pixels having intensity values greater than the first preliminary threshold value may be removed from further consideration. Pixels having intensity values less than the second preliminary threshold value may be removed from further consideration.

Referring to FIG. 4, a resultant image 405 (e.g., a first resultant image) may be formed by comparing image 410 (e.g., a first image) from a chemical band EM radiation detector (e.g., the first EM radiation detector) with image 420 (e.g., a second image) from a reference band EM radiation detector (e.g., the second EM radiation detector), the images having been subjected to the differential image quality enhancement methods described herein. In one or more embodiments, a plurality of resultant images (e.g., a first resultant image, a second resultant image, a third resultant image, etc.) may be analyzed. The plurality of resultant images may be formed from a plurality of images from the chemical band EM radiation detector (e.g., a first image, a third image, a fifth image, etc.) and a plurality of images from the reference band EM radiation detector (e.g., a second image, a fourth image, a sixth image, etc.). Determining the presence or absence of a chemical species may be based, at least in part, on one or more resultant images. A resultant image may be formed in any suitable manner, for example by determining the differences between intensity values of corresponding pixels of an image from a chemical band EM radiation detector and an image from a reference band EM radiation detector. One or more regions of interest may be identified based, at least in part, on non-zero intensity values for the resultant pixels. Differences observed in a resultant image may be considered a “positive” reading, i.e., an indication of the potential presence of a chemical species.

In one or more embodiments, the difference in intensity values of a resultant image may be formed by subtracting pixel intensities of an image from a reference band EM radiation detector (e.g., the second EM radiation detector) from pixel intensities of an image from a chemical band EM radiation detector (e.g., the first EM radiation detector). Subtracting pixel intensities of images yields a differential resultant image capturing a potential chemical species detected. One or more substantially spatially and/or temporally aligned pixels from the images (e.g., the first image and the second image) may be subtracted and the resulting intensity value used for the associated pixel in the resultant image (e.g., a first resultant image). The difference in intensity values in a resultant image may be represented by both the sign and the magnitude or may be absolute values.

In one or more other embodiments, the difference in intensity values of a resultant image may be determined by calculating the ratio between pixel intensities of an image from a reference band EM radiation detector (e.g., the second EM radiation detector) and pixel intensities of an image from a chemical band EM radiation detector (e.g., the first EM radiation detector). Calculating the ratio of pixel intensities of images yields a differential resultant image capturing any potential chemical species detected. The pixels of the images may be substantially spatially and/or temporally aligned. In one or more embodiments, the ratio is calculated between a plurality of pixels of an image from a chemical band EM radiation detector (e.g., the first EM radiation detector) and a plurality of pixels of an image from a reference band EM radiation detector (e.g., the second EM radiation detector). The ratio values are used to represent the intensity value for resultant pixels of the corresponding resultant image.

In one or more embodiments, a resultant threshold criteria may be applied to the resultant pixel intensity values. At least a portion of the resultant pixel values satisfying the resultant threshold criteria may be grouped into one or more regions of interest. The resultant threshold criteria includes at least one threshold value. In one or more embodiments, resultant pixels having an intensity value greater than a threshold value, and optionally additionally having the same sign, may be considered for inclusion in a region of interest. Taking into consideration the sign of the intensity of the resultant pixels may be used to remove resultant pixels from consideration due to low contrast, high noise such as those produced by vegetative backgrounds or highly reflective surfaces.

In one or more embodiments, the detection of a chemical species may be further enhanced using correlation coefficients, as described in U.S. Provisional Application Ser. No. 62/015,920, filed Jun. 23, 2014, and titled, “Method and Systems for Detecting a Chemical Species”, which is incorporated herein by reference in its entirety. For each of the one or more regions of interest identified in a resultant image, correlation coefficients may be determined. Each of the one or more regions of interest includes a subset of the plurality of pixels of the image. For example, one of the one or more regions may include a first subset of pixels from a plurality of pixels of an image from a chemical band EM radiation detector and a first subset of pixels from a plurality of pixels of an image from a reference band EM radiation detector and a first subset of pixels from a resultant image generated from the images. Another of the one or more regions of interest may include a second subset of pixels from the plurality of pixels of the image from the chemical band EM radiation detector and a second subset of pixels from the plurality of pixels of the image from the reference band EM radiation detector and a second subset of pixels from the resultant image generated from the same images. Additional regions of interest may include third subsets of pixels, fourth subsets of pixels, etc. Correlation coefficients (e.g., a first correlation coefficient, a second correlation coefficient, and a third correlation coefficient) may be calculated for each of the one or more regions of interest using the corresponding subset of pixels of the images. In one or more embodiments, the one or more regions of interest may be substantially the same as or different from the one or more areas of interest identified for any differential image quality enhancement methods described herein.

A correlation coefficient may be calculated using any suitable method. In one or more embodiments, the correlation coefficients may be calculated using a root mean square (RMS) correlation method to determine the correlation between corresponding pixel intensities within a region of interest to be correlated (e.g., the first image, the second image, and the first resultant image; the third image, the fourth image, and the second resultant image; the fifth image, the sixth image, and the third resultant image; etc.). A greater correlation between images of the chemical band EM radiation detector and the corresponding resultant image (e.g., first correlation coefficients) versus the correlation between images of the reference band EM radiation detector and the corresponding resultant image (e.g., second correlation coefficients) and optionally the correlation between images of the chemical band EM radiation detector and images of the reference band EM radiation detector (e.g., third correlation coefficients) may indicate the presence of a chemical species for a region of interest. In one or more embodiments, a plurality of regions of interest may be identified for analysis in a resultant image. Correlation coefficients for each of the plurality of regions of interest may be calculated using the intensity values of the corresponding pixels from an image from the chemical band EM radiation detector, an image from the reference band EM radiation detector, and the resultant image, similar to the calculations for the first, second, and third correlation coefficients of the first region of interest. Different subsets of the plurality of pixels corresponding to the regions of interest may be used to determine the correlation between images for the particular region of interest.

In one or more embodiments, the code of the analysis system may additionally be configured to direct the processor to: generate at least one resultant image, the resultant image including a plurality of resultant pixels, each pixel having an associated intensity value; determine one or more regions of interest based, at least in part, on the at least one resultant image; generate a plurality of correlation coefficients for at least one of the one or more regions of interest; and determine the presence of at least one chemical species based, at least in part, on the plurality of correlation coefficients. The detection method may additionally include generating a resultant image (e.g., at least a first resultant image); determining one or more regions of interest based, at least in part, on the resultant image; and determining correlation between an image (e.g., a first image) from the chemical band EM radiation detector (e.g., the first EM radiation detector), an image (e.g., the second image) from the reference band EM radiation detector (e.g., the second EM radiation detector), and the resultant image (e.g., the first resultant image) for at least one of the one or more regions of interest using a correlation coefficient algorithm to: calculate a first correlation coefficient, a second correlation coefficient, and optionally a third correlation coefficient.

The first correlation coefficient may be calculated using the intensity values of a first subset of pixels from the first plurality of pixels and a corresponding subset of resultant pixels from the plurality of resultant pixels within one of the one or more regions of interest. The second correlation coefficient may be calculated using the intensity values of a first subset of pixels from the second plurality of pixels and the corresponding subset of resultant pixels from the plurality of resultant pixels within the region of interest. The third correlation coefficient may be calculated using the intensity values of the first subset of pixels from the first plurality of pixels and the first subset of pixels from the second plurality of pixels within the corresponding region of interest. Additional correlation coefficients may be calculated for each of the additional regions of interest of the first image, the second image, and the first resultant image (e.g., a first correlation coefficient for an additional region of interest based on a second subset of pixels from the first image and the corresponding resultant pixels from the first resultant image; a second correlation coefficient for the additional region of interest based on a second subset of pixels from the second image and the corresponding resultant pixels from the first resultant image; and a third correlation coefficient for the additional region of interest based on the second subset of pixels from the first image and the second subset of pixels from the second image; etc.).

FIG. 5 illustrates flow chart for a method of detecting a chemical species using a detection system including multiple EM radiation detectors in accordance with one or more embodiments of the present disclosure. The detection method 500 begins at block 505. A beam of EM radiation is split into a plurality of beams (e.g., a first EM radiation beam and a second EM radiation beam) at block 510. At least a portion of one of the plurality of beams (e.g., the first EM radiation beam) is passed through a chemical band filter (e.g., a first bandpass filter) forming a filtered chemical band EM radiation beam (e.g., a filtered first EM radiation beam) at block 520 and at least a portion of another of the plurality of EM radiation beams (e.g., the second EM radiation beam) is passed through a reference bandpass filter (e.g., a second bandpass filter) forming a filtered reference band EM radiation beam (e.g., a filtered second EM radiation beam) at block 530. At least a portion of the filtered chemical band EM radiation beam (e.g., the filtered first EM radiation beam) is received by the chemical band EM radiation detector (e.g., the first EM radiation detector) at block 540 and at least a portion of the filtered reference band EM radiation beam (e.g., the filtered second EM radiation beam) is received by the reference band EM radiation detector (e.g., the second EM radiation detector) at block 550.

A differential image quality enhancement method is applied to images obtained from the chemical band EM radiation detector and the reference band EM radiation detector. In particular, at block 560, an image is obtained from the chemical band EM radiation detector (e.g., the first EM radiation detector). At block 565, an image is obtained from a reference band EM radiation detector (e.g., the second EM radiation detector). At block 570, one or more intensity values of a plurality of pixels of the image from the reference band EM radiation detector are adjusted based on one or more intensity value parameters from an image from the chemical band EM radiation detector according to one or more embodiments described herein. At block 575, one or more intensity values of a plurality of pixels of the image from the chemical band EM radiation detector are adjusted based on the one or more intensity value parameters from the image from the chemical band EM radiation detector according to one or more embodiments described herein. The adjustment of the intensity values of the reference band image and the chemical band image may be based on the images themselves or may be adjustment values determined from previous images from the chemical band and reference band EM radiation detectors (e.g., the first and second EM radiation detectors).

Once the differential image quality enhancement of the chemical band image and the reference band image is completed, at block 580 a resultant image (e.g., a first resultant image) is generated based on the chemical band image (e.g., a first image, a third image, a fifth image, etc.) and the reference band image (e.g., a second image, a fourth image, a sixth image, etc.). The resultant image is based, at least in part, on the data generated by the chemical band and reference band detectors (e.g., the first and second EM radiation detectors) according to one or more embodiments described herein. The detectors may be configured to generate a single image, a series of still images, or video images upon a command received from a computer system, a portion of which may be used to generate resultant images. At block 585, the presence or absence of a chemical species is determined based, at least in part, on one or more resultant images according to one or more embodiments described herein. If additional enhancement is desired, the method may return to block 510; otherwise, the enhancement process may end at block 590. Although not shown, a plurality of images may be generated and used to enhance the image quality, provide greater accuracy in the identification of a potential release, and decrease false positives resulting from background and noise interference.

A number of variations to the detection analysis techniques disclosed herein may be used to improve the reliability, ease of use, or ease of implementation of the detection system. In one or more embodiments, leak modeling results, leak detection criteria, camera and lens characteristics, and algorithm requirements may be combined to form deployment reference charts for setting up an autonomous detection system.

In one or more embodiments, additional detection analysis methods may be utilized to enhance the determination of the presence or absence of a chemical species and/or to increase the accuracy of the analysis. The additional detection analysis methods may be applied to chemical band images, reference band images, resultant images, and any combinations thereof. Additional analysis methods are described herein and further described in U.S. Patent Publication No. 2014/0002639 and U.S. Patent Publication No. 2014/0002667, which are incorporated herein by reference in their entirety.

Image analysis methods may use an algorithm to further analyze the images to distinguish chemical species (e.g., a gaseous plume of a chemical species) from other features in an image to improve accuracy and decrease the probability of false alarms. The algorithm may further distinguish the chemical species, such as hydrocarbon gas or vapor, from other ambient factors, such as water flows, steam plumes, furnace off gases, vehicles, persons, wildlife, and the like. In addition to using resultant images to determine the presence of a chemical species, enhanced identification techniques may be used. Such enhanced identification techniques may include analyzing for features such as deterministic features, probabilistic features, auxiliary features, and any combinations thereof. The image(s) may be in color or in grayscale, in which the difference in contrast may be used to identify features.

An additional detection analysis method may include analyzing for a deterministic feature. Deterministic features may include both spatial and kinematic features, among others. For example, the additional analysis method may determine geometric features, including the shape of a chemical plume and/or the size of a chemical plume, among others. The analysis may also determine shape constraints such as aspect ratio, dispersiveness (e.g., the thickness of the plume as a function of distance), convexity, and histogram of orientation gradient (HOG) of contour, among others. These features serve as constraints and provide screening of the potential regions which may represent a chemical species.

Kinematic or motion features may include determining that a plume is constantly moving but that the motion is restricted to a constrained area, as may be expected by a plume originating from a leak. Kinematic features may include size constraints of a plume, such as a minimal and maximal size through a sequence of images. The kinematic features may be used to filter out most rigid body interferences.

Probabilistic features may include a spatial pattern of the chemical plume, a temporal pattern of the chemical plume, or any number of other features. The additional analysis method may include joint spatial and temporal analyses such as a fast dynamic texture algorithm. In the probabilistic analysis, a statistical model described by two types of equations, e.g., evolution equations and observation equations, which respectively model the way the intrinsic state evolves with time and the way the intrinsic state projects to image pixels, may be fitted to the segmented pixel data. Parameters may be estimated by matrices. Other probabilistic analysis techniques may also be used, such as principal component analysis (PCA). In PCA, a determination of the variables causing changes to a plume is made, such as a statistical comparison of wind speed and direction with changes seen in plumes.

Auxiliary features may include comparing images visible to the human eye to the plume identified using the non-visible images, such as images in the IR spectrum. For example, the visible images may be used to differentiate organic vapor plumes and water steam. Generally, organic plumes may be dark in the non-visible images and not very visible in the visible images. In contrast, a steam plume may be bright in the non-visible images, due to emitted heat, and visible in the visible images. In addition to improving the detection, the visible images may be used to locate the leak in the plant environment, for example by comparing a registered image from a detector in the IR spectrum with an overlapping image from a detector in the visible spectrum.

The detection and confirmation of plumes may be enhanced by meteorological measurements collected by a meteorological monitor. The meteorological monitor may collect data on environmental conditions such as wind speed, temperature, precipitation, atmospheric haze, and the like. This data may then be used in one or more embodiments to confirm that a detected plume is consistent with the collected data. For example, the calculated motion of the plume may be compared with the wind direction, such as in a PCA algorithm. If the motion of the plume is inconsistent with the wind direction, the plume identification may be incorrect.

The code of the analysis system may be further configured to direct the processor to provide additional detection analysis, for example configured to direct the processor to analyze one or more images (e.g., chemical band images, reference band images, resultant images, and any combinations thereof) for deterministic features, probabilistic features, auxiliary features, and any combinations thereof.

In one or more embodiments, the detection reliability and/or accuracy may also be improved by utilizing chemical markers in various hydrocarbon streams. The chemical markers may be substances added to increase an absorbance or emission at a particular wavelength. Such markers may make the use of other detection techniques more effective. For example, fluorescent chemicals may be added to a hydrocarbon stream in very small amounts, such as a few parts-per-million, as these compounds often have a high quantum yield, which is the number of photons emitted divided by the number of photons absorbed. As the wavelength of light emitted may not overlap with natural sources, the identification of a plume from the fluorescence may be straightforward.

In one or more embodiments, EM radiation detectors having different bit depths and dynamic intensity ranges may be used. In such situations, the values applied to pixels of one of the EM radiation detectors may be adjusted proportionally.

It should be understood that the preceding is merely a detailed description of specific embodiments of the invention and that numerous changes, modifications, and alternatives to the disclosed embodiments can be made in accordance with the disclosure here without departing from the scope of the invention. Although the detection system may be described herein with reference to a first, chemical band EM radiation detector and a second, reference band EM radiation detector, one skilled in the art in light of the present disclosure will appreciate that any number of additional detectors may be used with the detection system described herein. The preceding description, therefore, is not meant to limit the scope of the invention. Rather, the scope of the invention is to be determined only by the appended claims and their equivalents. It is also contemplated that structures and features embodied in the present examples can be altered, rearranged, substituted, deleted, duplicated, combined, or added to each other. The articles “the”, “a” and “an” are not necessarily limited to mean only one, but rather are inclusive and open-ended so as to include, optionally, multiple such elements. 

What is claimed is:
 1. A method for differential image quality enhancement for a detection system including multiple electromagnetic radiation detectors, the method comprising: obtaining at least a first image from a first electromagnetic radiation detector, the image including a first plurality of pixels, each pixel having an associated intensity value; obtaining at least a second image from a second electromagnetic radiation detector, the image including a second plurality of pixels, each pixel having an associated intensity value; and adjusting one or more intensity values of the second plurality of pixels of the second image based on one or more intensity value parameters of the first image.
 2. The method of claim 1 further comprising adjusting one or more intensity values of the first plurality of pixels of the first image based on the one or more intensity value parameters for the first image.
 3. The method of claim 1, wherein adjusting the one or more intensity values of the second plurality of pixels includes adjusting the gain and offset.
 4. The method of claim 1, wherein the one or more intensity value parameters of the first image includes a maximum intensity (I_(max)) and a minimum intensity (I_(min)).
 5. The method of claim 1, wherein the adjusting of intensity values occurs periodically.
 6. The method of claim 1 further comprising: generating a histogram using at least a portion of the intensity values of the first plurality of pixels of the first image and at least a portion of the intensity values of the second plurality of pixels of the second image; adjusting the intensity values to offset either the second plurality of pixels to align with the intensity distribution of the first plurality of pixels or the first plurality of pixels to align with the intensity distribution of the second plurality of pixels, wherein the adjusting the one or more intensity values of the second plurality of pixels based on one or more intensity value parameters of the first image includes adjusting the gain of the second plurality of pixels and the one or more intensity value parameters of the first image includes a maximum intensity (I_(max)) and a minimum intensity (I_(min)).
 7. The method of claim 6 further comprising: identifying at least a first area of interest including one or more of the first plurality of pixels and the corresponding one or more of the second plurality of pixels that are substantially temporally and spatially aligned; and adjusting the intensity values to offset either the second plurality of pixels within the first area of interest to align with the intensity distribution of the first plurality of pixels within the first area of interest or the first plurality of pixels within the first area of interest to align with the intensity distribution of the second plurality of pixels within the first area of interest, wherein the adjusting the one or more intensity values of the second plurality of pixels based on one or more intensity value parameters of the first image includes adjusting the gain of the second plurality of pixels within the first area of interest and the one or more intensity value parameters of the first image includes a maximum intensity (I_(max)) and a minimum intensity (I_(min)) within the first area of interest.
 8. The method of claim 7, wherein the first area of interest consists of one or more pixels having intensity values within a defined range of intensity values.
 9. The method of claim 7, wherein the pixels within the first area of interest are determined by: identifying at least one of the first plurality of pixels and one or more adjacent pixels.
 10. The method of claim 7, further comprising: identifying a second area of interest including one or more other pixels of the first plurality of pixels and the corresponding one or more of the second plurality of pixels that are substantially temporally and spatially aligned; and adjusting the intensity values to offset either the second plurality of pixels within the second area of interest to align with the intensity distribution of the first plurality of pixels within the second area of interest or the first plurality of pixels within the second area of interest to align with the intensity distribution of the second plurality of pixels within the second area of interest, wherein the adjusting the one or more intensity values of the second plurality of pixels based on one or more intensity value parameters of the first image includes adjusting the gain of the second plurality of pixels within the second area of interest and the one or more intensity value parameters of the first image includes a maximum intensity (I_(max)) and a minimum intensity (I_(min)) within the second area of interest.
 11. The method claim 10, wherein the first area of interest and the second area of interest are at least a 10×10 grid of pixels.
 12. The method of claim 7, wherein the one or more of the first plurality of pixels within the first area of interest is determined by identifying a first pixel of the first plurality of pixels and one or more adjacent pixels.
 13. The method of claim 12, wherein the first pixel is identified by determining the differences between the intensity values of the first plurality of pixels and the corresponding second plurality of pixels and applying a threshold value to the determined differences in intensity values, the first pixel is selected from the differential intensity values greater than the threshold value.
 14. The method of claim 13, wherein the differences between intensity values of the first plurality of pixels and the corresponding second plurality of pixels is determined by calculating the ratio of the intensity values between each of the first plurality of pixels of the first image and the corresponding second plurality of pixels of the second image.
 15. The method of claim 13, wherein the differences between intensity values of the first plurality of pixels and the corresponding second plurality of pixels is determined by subtracting the intensity values between each of the first plurality of pixels of the first image and the corresponding second plurality of pixels of the second image.
 16. The method of claim 1 further comprising identifying one or more areas of interest, wherein the one or more areas of interest are determined by: obtaining a first baseline image from the first electromagnetic radiation detector, the first baseline image including a plurality of pixels, each pixel having an associated intensity value; obtaining a second baseline image from the second electromagnetic radiation detector, the second baseline image including a plurality of pixels, each pixel having an associated intensity value; determining the difference in intensity values between each of the plurality of pixels of the first baseline image and the plurality of pixels of the second baseline image, each of the plurality of pixels of the first baseline image and the plurality of pixels of the second baseline image are substantially temporally and spatially aligned; generating a first table of baseline information, the first table including the differences in intensity values between each of the substantially temporally and spatially aligned plurality of pixels of the first baseline image and the plurality of pixels of the second baseline image; determining the difference in intensity values between each of the first plurality of pixels of the first image and the second plurality of pixels of the second image, each of the first plurality of pixels and the second plurality of pixels are substantially temporally and spatially aligned; and applying the differences in intensity values of the first table to the difference in intensity values between each of the substantially temporally and spatially aligned first plurality of pixels and second plurality of pixels, the pixels proximate one another representing non-zero values are grouped to form the one or more areas of interest, wherein the adjusting the intensity values of the second plurality of pixels occurs within each of the one or more areas of interest based on one or more intensity value parameters of the first image for each of the one or more areas of interest.
 17. The method of claim 16, wherein the one or more intensity value parameters of the first image includes a maximum intensity (I_(max)) and a minimum intensity (I_(min)).
 18. The method of claim 16 further comprising: obtaining a third baseline image from the first electromagnetic radiation detector, the third baseline image including a plurality of pixels, each pixel having an associated intensity value; obtaining a fourth image from the second electromagnetic radiation detector, the fourth image including a plurality of pixels, each pixel having an associated intensity value; determining the difference in intensity values between each of the plurality of pixels of the third baseline image and the plurality of pixels of the fourth baseline image, each of the plurality of pixels of the third baseline image and the plurality of pixels of the fourth baseline image are substantially temporally and spatially aligned; determining a moving average of the differences in intensity values of the first baseline image, the second baseline image, the third baseline image, and the fourth baseline image for use as the baseline information; determining the difference in intensity values between each of the first plurality of pixels of the first image and the second plurality of pixels of the second image, each of the first plurality of pixels and the second plurality of pixels are substantially temporally and spatially aligned; applying the differences in intensity values of the moving average baseline information to the difference in intensity values between each of the first plurality of pixels and the second plurality of pixels, the pixels proximate one another representing non-zero values are grouped to form one or more areas of interest; and adjusting one or more intensity values of the second plurality of pixels within each of the one or more areas of interest based on one or more intensity value parameters of the first image for each of the one or more areas of interest.
 19. The method of claim 16, wherein the difference in intensity values is calculated by subtracting the intensity values.
 20. The method of claim 16, wherein the difference in intensity values is determined by calculating the ratio of the intensity values.
 21. A system for detecting at least one chemical species comprising: a lens; a beam splitter; a first electromagnetic radiation detector; a second electromagnetic radiation detector; and an analysis system including code within a processor, a non-transitory, computer-readable medium, or a combination thereof, the code configured to direct the processor to: identify at least a first image from the first electromagnetic radiation detector, the first image including a first plurality of pixels, each pixel having an associated intensity value, identify at least a second image from the second electromagnetic radiation detector, the second image including a second plurality of pixels, each pixel having an associated intensity value, and adjust one or more intensity values of the second plurality of pixels of the second image based on one or more intensity value parameters of the first image.
 22. The system of claim 21, wherein the analysis system automatically adjusts the one or more intensity values.
 23. A method of detecting at least one chemical species comprising: providing a detection system including a first electromagnetic radiation detector and a second electromagnetic radiation detector; splitting an electromagnetic radiation beam into at least a first electromagnetic radiation beam and a second electromagnetic radiation beam; receiving at least a portion of the first electromagnetic radiation beam by the first electromagnetic radiation detector; receiving at least a portion of the second electromagnetic radiation beam by the second electromagnetic radiation detector; generating at least a first image by the first electromagnetic radiation detector and at least a second image by the second electromagnetic radiation detector; enhancing the image quality of the first image and the second image according to the method of claim 1; and determining the presence of at least one chemical species based, at least in part, on a comparison of the enhanced first image and the enhanced second image. 